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Application of Darboux Transformation to solve Multisoliton Solution on Non-linear 

Schrodinger Equation 
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Darboux transformation is one of the methods used in solving nonlinear evolution equation. Ba- 
sically, the Darboux transformation is a linear algebra formulation of the solutions of the Zakharov- 
Shabat system of equations associated with the nonlinear evolution equation. In this work, the 
evolution of monochromatic electromagnetic wave in a nonlinear-dispersive optical medium is con- 
sidered. Using the Darboux transformation, explicit multisoliton solutions (one to three soliton 
solutions) are obtained from a trivial initial solution. 

PACS numbers: 



I. INTRODUCTION 

In 1882, Gaston Darboux introduced a method to solve 
Sturm-Liouvillc differential equation, which is called Dar- 
boux transformation afterwards [lj . Remarkably, this 
transformation can solve not only Sturm-Liouville dif- 
ferential equation (another form of time-independent 
Schrodinger equation), but also for the other forms 
of linear and non-linear differential equations, such 
as Korteweig-de Vries, Kadomtsev-Petviashvilli, Sine- 
Gordon, and Non-linear Schrodinger equations. Fol- 
lowing exposition explains about Darboux transforma- 
tion and its generalized form called Crum's theorem on 
Sturm-Liouvillc differential equation. 



II. DARBOUX TRANSFORMATION AND 
CRUM THEOREM 

Consider following Sturm-Liouvillc differential equa- 
tion 



Ipxx + Wp = Xljj 



(1) 



where u is function of x and A is a constant. In the 
Schrodinger equation, u represent a potential. Darboux 
transformation is defined as follow Q 



<A[1] = (J^ - <rij ^ = il>x - ^V> 



W (ip,ip\l\) expresses Wronskian determinant as below 



(2) 



W(yji,tpi, ...il) N ) = 


4>i 


4>2 


IpN 

■ 


■ (3) 




4 n ~ 1] 


4 nil) 


"■■ 4r 1} 





where tp^ means ip derivative for rt-times and ipi is the 
solution of ip for A = Ai. If ip is a solution, ip[l) is the so- 



lution for the following Sturm-Liouville differential equa- 
tion then 



ij xx [l}+u[l}ip[l} = \ip[l} 



(4) 



where u[l] is a new function of u which has been trans- 
formed. It will be shown below that Darboux transfor- 
mation acting on ^[1] influences the potential u, if ip[l] 
is invariant over Eq. ©. 

From Eq.©, one can obtain following relation 

- rp xx [l] = -tpxxx + o lxx ip + 2<r lx ip x + o\ip xx (5) 

Substitution of Eq.© and Eq.© into Eq.© results 

- ip xxx + <J lxx ip + 2<r lx ip x + aii> xx + u[l] (ip x - axip) 

= A (ip x — a-iip){6) 

Using Eq.© to substitute ip xx results 



(u[l] -u + 2ai x )ip x + 

x + <J\U — <T\U\\\) Ip = 



(7) 



It is clear from Eq.©, the following relations can be 
obtained 



u[l] = u - 2o\ x 
o-ixx - u x + 2u\U\ x = 



(8) 
(9) 



The Eq. © shows that function u is also under the trans- 
formation due to ip in Eq. © under Darboux transforma- 
tion. In other words, Sturm-Liouville equation in Eq.© 
is covariant under Darboux transformation action 

ip —> ipm an d u — ► u[l] 

Interestingly, Darboux transformation can be recur- 
sively applied to Sturm-Liouville equation solution and 
consequently the potential is under transformation to en- 
sure that the solution belongs to the equation. Below will 
be discussed the consequence if the solution in Eq. © is 
acted by Darboux transformation twice. 



, d lp2x[l] 
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where 



Therefore 



0[2] = lp2x ~ -^02 
01 



(11) 
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01 (ipxxfax - 1p2xx1px) - 02 {^xx^lx - 4>lxx4>x) + ~0 {^2xxj>\x ~ 01aa02z) 



01 "02 
01x ^22; 0a 

XX 



W(0i,0 2 ) 



(12) 



From the definition i[>[l] and 0[2], u[l], and u[2], Crum 
proposed general form of Darboux transformation in the 
following form 



W(ip 1: ...,ip N ,ip) 



(13a) 



'-)- 

u[N}=u-2— lnW{iP u ...,ip N ) (13b) 



where tp[N] satisfies the equation 

-iP xx [N]+u[N]i)[N] = \^[N] 



(14) 



The Crum's theorem is based on the fact that Darboux 
transformation acting iV-times on a function can be writ- 
ten 

ip[N] = D[N]ip = ip (N) + Sl ^ N -V + ... + SN ip. (15) 

Ansatz of Crum's theorem is that sn is obtained from 
the following relation 



N 



N-k 
j 



-0 



(N) 



(16) 



fe=i 



where j=l, 2, ... , N. The equation (UHJ) can be explained 
and arranged in matrix representation by following 



/ 4 N -^ 




( \ 



(17) 



By the use of Cramer's rule, one can obtain Sj for j = 1, 



, N. For instance, si 



si 



-■01 ¥>1 01 



< 0^ 



(JV-1) 



0i 







01 



For example, the implementation of Crum's theorem as 
given above will be proven below for N=2. The equation 
(fI7J)) for iV=2 can be explained by following 



siipix + s 2 ipi = -Iplxx 

Sli>2x + S 2 1p2 = ~4>2xx 



(18a) 
(18b) 



By the use of Cramer's rule, one can determine s\ and 
S2 as below 



-'ipixx 0i 




01 


Iplxx 


~lp2xx 02 




"02 


1p2xx 


Iplx Ipl 




01 


02 


1p2x 02 




Iplx 


1p2x 



(19a) 



si 



Iplx - 


■Iplxx 




010! 


^IXX 


4>2x ~ 


■1p2xx 




02z 


4>2xx 




Ax 


01 






01 


02 




1p2x 


02 






Iplx 


02z 



(19b) 



Darboux operator in Eq. (fl5|) for N = 2 can be described 
0[2] = D[2]0 = tp xx + si0 x + s 2 0. (20) 
Whilst, substitution Eq.([T9]) into Eq.® yields 

1^(01,02,0) 



0[2] = £[2]0 = 



^(0i,0 2 ) 
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III. DARBOUX TRANSFORMATION ON 
NONLINEAR SCHRODINGER EQUATION 

From Zakharov-Shabat system, one can derive sev- 
eral nonlinear differential equations, such as nonlinear 
Schrodinger and Sine-Gordon equations [l[. In 1982, 
Bobenko and Salle can apply Darboux transformation 
on generalized Zakahrov-Shabat system 0. In this sec- 
tion, the Darboux transformation application for first or- 
der linear Zakharov-Shabat system will be explained, fol- 
lowed by its application for nonlinear Schrodinger equa- 
tion. 



A. Darboux Transformation for First Order Linear 
System 

Consider the following system: 

ip x = —iXip + iqcj>. (21a) 

4> x = -X<f> + inp. (21b) 

These differential equations can be assembled into the 
following matrix: 



(22) 



where 

J = 



-i 
i 



A 



Ai 
A 2 



U = 



iq 
ir 

4>\ 4>2 



The solution for Ea.([2"T]) is 



4>k 

<t>k 



where k=l, 2. Dar- 



boux transformation acting on Eo. ([2"2"|) is defined by fol- 
lowing 

* -> = D[l]# = *A + Si^ = *A - (^ 1 A 1 ^ 1 ) 

(23) 

where is the solution obtained for Eq. (|2"2"|) with A = 
Ai. The equation (|2"2")l is called covariant under Darboux 
transformation if the potential is transformed 



U[1] = U + [S 1 ,J]. 



(24) 



Therefore, is the solution of following differential 

equation 



(25) 



*[1] = J*[1]A + [/[!]*[!]. 



Acting iV-timcs Darboux transformation on Eq. (|22|) . it 
means following this rule 



D[N]V = *[JV] = *A iV +S , i1'A 



,JV-1 



-Sn-^A+Sn^- 
(26) 



Whilst, generalized form of the potential transforms 

U[N] = U+[S N ,J}. (27) 

From Crum's theorem, Si can be obtained from the fol- 
lowing condition 

*[iV] A=Afc , *=* fc = 0. (28) 

Because k is running until N, the set of equations in 
Eg. ([28]) can be rewritten in the following form then 



,N-1 



Si*ia; v 
I 

S^nA?- 1 - 



,N-2 



I 



\-S N *i = -*iAf . 

I i 

-S N V N =-*jvAf. 

(29) 

If the equations in (|29p are assembled into matrix form 



5 a *iAi v 
I 

S 2 V N A?- 2 - 



Si 



Sn 



VnA^ 1 



( -*iA{ 



N 



Where 

Sn ■ 



Cl2 
N 
C22 
°N 



V N = 

A N 



-VnA n , . 

1p2N-l ^p2N 
02JV-1 4>2N 

A27V-1 
X2N 



(30) 



Si, for i=l, ... , N, can be determined by the use of 
Cramer's rule. Therefore, for Si, one can obtain 



si 1 si 2 



(31) 



Sf S 2 

-V^Af-V^Af" 2 ^,...,^,^), 



where 

A = det(Af" 
S? 

Sf = det{K 
Si 2 = det(Af-V i ,Af^,Af- 2 ^,Af- 2 (/) ! :,...,^,^), 
S 21 = det(A.f^,Af- 1 (/) ? :,Af- 2 V ! :,Af" 2 ^,-,^,^)- 
Where i=l, ... , 2N and also using following definition 



det(AfVi)Af Vij^f 2 V'i,Af 2 c/)i,...,ipi,<j)i), 

'^xf^x^-^xf- 2 ^,...,^,^), 





/ ai 61 . 




det (a, ; , 6j, /,;) = det 




• h 




c ::: : 


: ::.) 



Consequently, iV-timcs Darboux transformation on 
Eq. (f2"T|) , q[N] and r[N] are expressed 



012 

q[N] = q + 2^ 



r[N] 



21 



' A 



(32a) 
(32b) 



The application of the entire calculation in this section 
for obtaining multisoliton of nonlinear Schrodinger equa- 
tion will be explained in the following section. 
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B. Darboux Transformation and Cram's Theorem 
on Nonlinear Schrodinger Equation 

Below, the derivation of nonlinear Schrodinger equa- 
tion from generalized linear Zakharov-Shabat equation 
will be given, followed by the explanation of Darboux 
operator implementation on its solution then. 

Consider following two linear equations 



^kx = J^kAk 



(33a) 



* fct = 2J«Z k A 2 + 2UV k A k + (JU 2 - JU X ) **. (33b) 



Where 



i>2k-i i>2k 



92fc-l <P2k 

From the integrability condition, ^ xt = ^tx, it yields 



(34) 



(J 2 U 2 - J 2 U X - 2U X + JU X J - JU 2 J) *A + 
(U f + UJU 2 - UJU X - 2JUU X - JU 2 U + JU XX + JU X U) * = 0. (35) 

I 



If Eq. ([35]) is assembled in matrix form 



qr - q x r x 
qr - q x r x 



*A 



-iqr x + 2iqr x - iqr x 
ir t + qr 2 + q x r x r + r xx 

I 



9 

q r 



iqr x 



qqxT x q xx I ^ q 

2iqr x + iqr x 



(36) 



For q = r* and r = q* , this consistency equation (|36p 
gives 

ir t + 2r*r 2 + r xx = 0, (37a) 

iq t - 2q*q 2 - q xx = 0. (37b) 
It is clear that the equation (|37p yields nonlinear 



Schrodinger equation. Furthermore, it will be shown the 
consequence if the two linear equations in (|33[) are trans- 
formed under Darboux transformation (|!i3"|) . ^[l] = 
J*[l]A+t/[l]*[l] or (*A - o^) x = J[l] (*A - Si*) A+ 
U[l] (*A - Sj*). The substitution of the equation l|3"3"]) 
results 



(J - J[l]) *A 2 + ([/■ - Si J + J[l]Si - £/[!]) *A + {S lx + SiU - U[l]St) * = (38) 



From the equation (|38|) of each coefficient, one can de- accomplished, 
termine 

C. Darboux Transformation Solution of Nonlinear 
J = J[l] (39a) Schrodinger Equation for Trivial Initial Solution 



U[l] = U + J[l]S 1 -SiJ =S> U[l] =U +[J, Si] (39b) 

S lx + U - U[1]S! = (39c) 

where J[l] and U[l] are respectively J and U in the 
equation (|33[) after Darboux transformation has been 



In this section, we will examine the consequence on 
nonlinear Schrodinger equation if its initial solution is 
trivial. For the case r = q = 0, the equation (|3"3"|) will be 
in the form 
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l/>2k-l lp2k 
4>2k-l 02k 



i W 1p2k-l *l>2k 
-1 I \ 4>2k-l <p2k 



A 2fe -i 
A 2fe 



(40a) 



i>2k-i *l>2k 

<p2k-l <hk 



= 2 



i 
-i 



1p2k-l i>2k 
4>2k-l <j>2k 



r 



A 2fe -i 
A 2fc 



(40b) 



The equations in (|40|) yield the following solution 

fok-i - A 2fe _ ie J ( A -- ia;+2A --0; ^ 2fe = A 2fc e< A2fc *+ 2A -) 

(41a) 

<hk-i = B^-ie-'t^-^+^-O; V. 2fc = S 2fce -< A — +2A -) 

(41b) 

The solutions in (j4"Tj) show that from the trivial solution 
in (|33p by using Darboux transformation can create mul- 
tisoliton as mentioned in the following sections. 



Furthermore, by choosing tp2 — — 0*, it is resulted 



A 



i>\ i>2 

4>1 02 



(^102 ~ 1p2^l) = (V'lV'l ~ 01 



Ivll 2 + I01l 2 



(45) 



IV. MULTISOLITON OF NONLINEAR 
SCHRODINGER EQUATION USING DARBOUX 
TRANSFORMATION 

In this chapter, the calculation of multisolition from 
trivial solution using Darboux transformation will be ex- 
plained. In the previous chapter has been shown that 
by choosing appropriate linear Zakharov-Shabat system 
and also from consistency condition, two types of nonlin- 
ear Schrodinger equation can be obtained, where either 
of them is 



in + r xx + 2 \r\ r = 0. 



(42) 



For iV-times Darboux transformations, the equation is 
matrix element of E/[./V], i.e. (U[N}) 21 , which is in the 
following form 



i(r[N)) t + (r) xx + 2\r[N]\ 2 r[N}=0. 



(43) 



Below, we will describe how to generate multisoliton so- 
lution from trivial solution for 3-times Darboux transfor- 
mation action. 



A. The Calculation of 1-Soliton of Nonlinear 
Schrodinger Equation using Darboux 
Transfer mat ion 



For the case N = 1, according to the previous expla- 
nation, from the equation (|3ip yields 



Sf 



Ai0i A 2 2 
4>i 4>2 



(Ai - A 2 ) 0i0 2 = (Ai - A*) i>X§\. 

(44) 



Consequently, the initial trivial solution of equation 
gives 



r[l] 



S? (Ai - \l)tpl4>i 

~ Z A '-<■ !2 , ,2 



(46) 



Let us define the following variable 

6 = Aia; + 2U 2 = Re (0) + ilm (0) . (47) 



Substitution the equation (|47[) into (|4Tj) for the case k = 
1, followed by substitution into (|4"6"| afterwards. One can 
determine 

(2iIm(A 1 ))A 1 B 1 e(- Mm ( e )) , . 

r 1 = v in ^ (48) 

(lA.l 2 e 2(Ro(f)) + l^l'e- 2 ^))] 

Furhtermore, let us define the following constants 



A 1 



(49a) 



(49b) 



Substitution the equation (j49f into (|47|) . one can obtain 
the following solution 

U cosh(2(Rc(0) + $ 2 )) V ; 

The equation (JSTJJ) is the solution of 1-soliton of nonlin- 
ear Schrodinger equation which has the form (|42p . From 
(|50p . it can be shown that imaginary component of Ai 
related to wave amplitude, while its real component rep- 
resents the wave velocity. 



B. The Calculation of 2-Soliton of Nonlinear 
Schrodinger Equation using Darboux 
Transformation 

For 2-soliton, according to the Crum's theorem 

Q21 



12] 



A 



(51) 



6 



Where 



21 



\ 2 fa \\fa 2 A§0 3 Ai</> 4 
Ai^i A 2 A 3 3 A 4 4 

^1 02 03 04 



(52a) 



((Re (Ai) - Re (A3))) - 1 (Im (Ai) - Im (A3)) = -p 

(54b) 



(Re (Ai) - Re (A 3 )) + i (Im (Ai) + Im (A 3 )) = (3 (54c) 



A 



Al01 A 2 ^ 2 A 3 3 A 4 4 

Al0l A2 A303 A404 

01 02 03 04 

01 fa fa fa 



(52b) 



(Re(Ai)-Re(A 3 ))-j(Im(Ai)+Im(A 3 )) =0*. (54d) 



Furthermore, by substituting ((52)) into (|5Tj) . one can ob- 
tain 

A?0i Al^>| \\fa Ai0 4 
Ai0i A 2 A 3 3 A 4 4 

01 02 03 04 
01 02 03 04 



12] 



Where 



6 = Aix + 2iA 2 = Re (0) + ilm (6>) (54e) 



(53) 



Al0l A 2 0>2 A 3 "03 A 4 04 

Ai0i A 2 A 3 3 A 4 4 

01 02 03 04 

01 02 03 04 



By defining the following parameters 
(Re (Ai) - Re (A3)) + i (Im (Ai) - Im (A 3 )) = p* (54a) 



a = A3X + 2t\ 2 3 = Re (a) + ilm (a) (54f) 



One can obtain the solution for 2-soliton as stated below 



r[2) = Ai 



Im(Ai)e Mni W (p*/3e- 2R °( Q ) +p/3*e- 2Rc < Q )) -Im(A 3 )e Mm ( Q ) (p* (3* e~ m °^ + pf3e- 2Rc ^) 



\p\ z (cosh (2 (Rc (6>) + Re (a)))) 



(cosh (2 (Re (6<) - Re (a)))) 



41m (Ai) Im (A 3 ) (cosh (2 (Im (6) + Im (a)))) 

(55) 



The equation (|55| is obtained by choosing fa = — 0>* and 

04 = -03- 



The Calculation of 3-Soliton of Nonlinear 
Schrodinger Equation using Darboux 
Transformation 



3-soliton can be generated by similar manner as ob- 
taining 2-soliton. First, determinant form of Sf 1 and A 
which are in the following forms 



A 



A?0i 
A 2 0i 
Ai0i 
Ai0i 
0i 



Ai0 2 

A|02 
A 2 02 
A 2 02 

02 

02 



Al.03 

a§03 

A 3 03 
A3 03 
03 



A^04 
A 4 04 
A 4 4 
A 4 4 

04 

04 



Af,05 

Ai0 5 

A 5 05 
A505 

05 

fa 



A 2 . 06 

Ai0 6 

A 6 06 
A606 

06 

06 



(56b) 



In addition, substitute (|56p to the following equation 

C21 



r[3] 



A 



(57) 



Xffa A|0 2 A|0 3 Al0 4 Ai0 5 A|0 6 



A?0i 
Ai0i 
Ai0i 
0i 

01 



A 2 02 
A 2 02 
02 

q> 2 



a§03 

A3 03 
A3 03 

03 

03 



A 4 04 
A 4 04 
A 4 04 

04 

<?>4 



Ai0 5 

A505 
A505 

05 

05 



A|06 

Ae06 

A606 
06 
00 



(56a) 



The equation (|57j) yields the solution for 3-soliton by 
obeying the following chosen 



02 = -(, 



"03 



(58a) 



(58b) 




50 25 -25 -50 



x 



(b)Contour. 
FIG. 1: 2 solitons. 



A = ~<t>l (58c) 

Below the illustration based on the above formulation 
for three solitons are given. By changing the parameters, 
several forms of multi-solitons can be obtained. 



D. Conclusion 

It has been shown that nonlinear refractive index in 
fiber optics materials, which has Kerr-like medium char- 
acteristics, will compensate dispersive effect of electro- 
magnetic wave propagating in that fiber optics materi- 
als. The nonlinear evolution equation of electromagnetic 
wave due to the both effects can be solved by Darboux 
transformations. The advantage of Darboux transforma- 
tion is that one can obtain multisolitons solution from 
trivial solution. This transformation has the similar form 
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(c)Contour of 3 Solitons. 
FIG. 2: 3 soliton contour 



with the other methods which is used to solve nonlinear 
differential equations, i.e. Cole-Hopf transformations. 
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